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Abstract. The three-dimensional equations of compressible magnetohydrodynamic 
isentropic flows are considered. An initial-boundary value problem is studied in a 
bounded domain with large data. The existence and large-time behavior of global weak 
solutions are established through a three-level approximation, energy estimates, and 
weak convergence for the adiabatic exponent 7 > § and constant viscosity coefficients. 



1. Introduction 

Magnetohydrodynamics (MHD) concerns the motion of conducting fluids in an electro- 
magnetic field with a very broad range of applications. The dynamic motion of the fluid and 
the magnetic field interact strongly on each other. The hydrodynamic and clectrodynamic 
effects are coupled. The equations of three-dimensional compressible magnetohydrody- 
namic flows in the isentropic case have the following form ( [2| \TE\ fT9] ) : 

Pt + div(pu) ^ 0, 

{pu)t + div (pu ® u) + Vp = (V X H) X H + pAu + (A + ^)V(divu), (1.1) 
Ht - V X (u X H) = -V X (j/V X H), divH = 0, 

where p denotes the density, u e M'^ the velocity, H e M"^ the magnetic field, p{p) = ap'^ the 
pressure with constant a > and the adiabatic exponent 7 > 1; the viscosity coefficients 
of the flow satisfy 2/i + 3A > and p > 0; ly > is the magnetic diffusivity acting as a 
magnetic diffusion coefficient of the magnetic field, and all these kinetic coefScients and the 
magnetic diffusivity are independent of the magnitude and direction of the magnetic field. 
The symbol (8) denotes the Kronecker tensor product. Usually, we refer to the first equation 
in (jl.ip as the continuity equation, and the second equation as the momentum balance 
equation. It is well-known that the electromagnetic fields are governed by the Maxwell's 
equations. In magnetohydrodynamics, the displacement current can be neglected f |18[|19j ). 
As a consequence, the last equation in p.l|) is called the induction equation, and the electric 
field can be written in terms of the magnetic field H and the velocity u, 

E = j/VxH-uxH. 

Although the electric field E does not appear in (jl.ip . it is indeed induced according to 
the above relation by the moving conductive flow in the magnetic field. 

In this paper, we are interested in the global existence and large-time behavior of so- 
lutions to the three-dimensional MHD equations (jl.ip in a bounded domain 51 C K."^ with 
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the following initial-boundary conditions: 



' p{x,Q)^ Po{x)eL'<{n), po(2:)>o, 

0)u(a::, 0) = mo(a;) € ^^(f]), mo = if po = 0, 
H(a;,0) = Ho(x) e L'^{il), divHo = in V'{^), 
u|an-0, H|aQ = 0. 




Gil 



(1.2) 



There have been a lot of studies on MHD by physicists and mathematicians because of 
its physical importance, complexity, rich phenomena, and mathematical challenges; see 
[31 El ISl 13 [111 [13 lining and the references cited therein. In particular, the one-dimensional 
problem has been studied in many papers, for examples, [21 SI [3 [13 [13 [HI US] and so 
on. However, many fundamental problems for MHD are still open. For example, even 
for the one-dimensional case, the global existence of classical solution to the full perfect 
MHD equations with large data remains unsolved when all the viscosity, heat conductiv- 
ity, and diffusivity coefficients are constant, although the corresponding problem for the 
Navier-Stokes equations was solved in [16] long time ago. The reason is that the presence 
of the magnetic field and its interaction with the hydrodynamic motion in the MHD flow 
of large oscillation cause serious difficulties. In this paper we consider the global weak 
solution to the three-dimensional MHD problem with large data, and investigate the fun- 
damental problems such as global existence and large-time behavior. A multi-dimensional 
nonisentropic MHD system for gaseous stars coupled with the Poisson equation is studied 
in |S], where all the viscosity coefficients depend on temperature, and the pressure depends 
on density asymptotically like the isentropic case p{p) ~ ap'^ . In this paper, we study 
the multi-dimensional isentropic problem p.ip - (|1.2p with 7 > |, where all the viscosity 
coefficients p, A, v are constant. We remark that 7 = | for the monoatomic gases. 

When there is no electromagnetic field, system (|l.ip reduces to the compressible Navier- 
Stokes equations. See [HllIllHI] and their references for the studies on the multi-dimensional 
Navier-Stokes equations. In particular, to overcome the difficulties of large oscillations of 
solutions, especially of density, the concept of a renormalized solutions is used in [211 [S] ■ 
Based on this idea, we study the initial-boundary value problem (|l.ip - (|1.2p for the MHD 
system in a bounded three-dimensional domain 51. The goal of this paper is to establish 
the existence of global weak solutions for large initial data in certain functional spaces 
for 7 > I and to study the large-time behavior of global weak solutions when the mag- 
netic field and interaction present. The existence of global weak solutions is proved by 
using the Faedo-Galerkin method and the vanishing viscosity method. We first obtain a 
priori estimates directly from ()1.1|) . which is the backbone of our result. In the proof of 
the existence, we use the similar approximation scheme to that in |10j which consists of 
Faedo-Galerkin approximation, artificial viscosity, and artificial pressure. Then, motivated 
by the work in [5] , we show that an improvement on the integrability of density can ensure 
the effectiveness and convergence of our approximation scheme. More specifically, we show 
that the uniform bound of p'^ln(l -|- p) in L^, rather than the uniform bound of p'^^^ in 
for some > as used in [lOl [9] [21], ensures the vanishing of artificial pressure and 
the strong convergence of the density. To overcome the difficulty arising from the possible 
large oscillations of the density p, we adopt the method in Lions [3T] and Feireisl [3] which 
is based on the celebrated weak continuity of the effective viscous flux p—{\ + 2/x)divu (see 
also Hoff [H]). The estimates obtained by our approach produce further the large-time 
behavior of the global weak solutions to the initial-boundary value problem (|l.ip - (|1.2p . To 
achieve our goal for the MHD problem, we also need to develop estimates to deal with 
the magnetic field and its coupling and interaction with the fluid variables. The nonlinear 
term (V x H) x H will be dealt with by the idea arising in incompressible Navier-Stokes 
equations. 
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We organize the rest of the paper as follows. In Section 2, we derive a priori estimates 
from (jl.ip . give the definition of the weak solutions, and also state our main results and 
give our approximation scheme. In Section 3, we will show the unique solvability of the 
magnetic field in terms of the velocity field. In Section 4, following the method in [9] and 
the result obtained in Section 3, we show the existence of solutions to the approximation 
system. In Section 5, we follow the technique in [lOj with some modifications to get the 
strong convergence of p in i^((0, T) x 51). In Section 6, motivated by [TT] we study the 
large-time behavior of global weak solutions to (|l.ip - (|1.2p . 



2. Main Results 

In this section, we reformulate the initial-boundary value problem (jl.ip - (jl.2[) and state 
the main results. 

We first formally derive the energy equation and some a priori estimates. Multiplying 
the second equation in (jl.ip by u, integrating over fi, and using the boundary condition 
in (11.21). we obtain 



(2.1) 



j j (\pn' + ^ Adx+ I (A.|i?u|2 + (A + /i)(divu)2)dx 

^ [ ((V X H) X H) • u dx. 

The term on the right hand side of (|2.1[) can be rewritten as 

^ ((V X H) X H) • u dx = - ^H^VuH + ^VdHp) • dx. 

Hence, (|2.ip becomes 

I / flpu^ + ^pAdx+ I (Ai|2?up + (A + M)(divu)2)dx 



H^VuH+ ivdHp) . u I dx. 



(2.2) 



Multiplying the third equation in (|l.ip by H, integrating over fJ, and using the boundary 
condition in (11.21) and the condition divH = 0, one has 



ld_ 

Direct calculations show that 



[ IHpdx + [ (V X (j/V X H)) • H dx = [ (V X (u X H)) • H dx. (2.3) 
Jn Jn Jn 

ow that 

/ (V X (i/V X H)) • H dx = / |V X Hpdx, 
Jn Jn 

^ (V X (u X H)) • H dx = ^ j^H^VuH ^V(|Hp) • dx. 



Thus dM]) yields 

1 d 
2dt 



[ IHpdx + z/ / |VxHpdx=: / (H^VuH + iv(|H|2).u)dx. 
Jn Jn Jn \ 2 / 



Adding (gj]) and ([23) gives 



/ (/i|-Du|2 + (A -I- n){divuf + X Hp) dx = 0. 
Jn 



(2.4) 



(2.5) 
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From Lemma 3.3 in |23j . our assmirptions on initial data, and (|2.5p . we have the following 
a priori estimates: 

vi(.L\[Q,T]-Hl{n))- 
■ReL\[0,T]-Hlmf^L^{[Q,T]-L\^))- 

peL^{[o,nL-'{n)). 

Multiplying the continuity equation (i.e., the first equation in p.ip ) by b' {p), we obtain 
the renormalized continuity equation: 

b{p)t + div(6(p)u) + {b'{p)p - 6(p))divu = 0, (2.6) 

for some suitable function b e C^(M+). Following the strategy in (51] we introduce 
the concept of finite energy weak solution (p, u, H) to the initial-boundary value problem 
(|l.ip - (|1.2p in the following sense: 

• The density p is a non-negative function, 

p e C{[0, T];L\n)) n L°°([0, T];L^n)), p{x, 0) po, 
and the momentum pu satisfies 

pueC([0,T];i£,(f})); 

• The velocity u and the magnetic field H satisfy the following: 

u&L\[0,T];H^in)), HGL2([0,T];i/i(r!))nC([0,r];iLafc(f^)), 
pu (g) u, V X (u X H), and (V x H) x H are integrablc on (0, T) x fl, and 
pu{x, 0) = mo, H(a;, 0) = Ho, divH = in V'{n); 

• The system (jl.ip is satisfied in P'(M'' x (0, T)) provided that p, u, and H are 
prolonged to be zero outside ft; 

• The continuity equation in is satisfied in the sense of renormalized solutions, 
that is, jlH) holds in V'{n x (0,T)) for any b G Ci(R+) satisfying 

6' (2;) = for all z e large enough, say, z > zq, (2.7) 

where the constant zq depends on the choice of function b; 

• The energy inequality 

I I (Ai|£'u|2 + (A + p)(divu)2 + jy|VxH|2)da;ds<i;(0), 
holds for a.e t E [0,T], where 

i.(.) = /^(i,u^ + -^p^ + i|Hp)d., 

and 

Eio)^ [ (l^^ + ^p2 + l\nAdx. 

Jn\2 Po 7-1 2 J 

Remark 2.1. As a matter of fact, the function b docs not need to be bounded. By Lebesgue 
Dominated convergence theorem, we can show that if p, u is a pair of finite energy weak 
solutions in the renormalized sense, they also satisfy (|2.6p for any b e C^(0, 00) n C[0, 00) 
satisfying 

|^'('Z)-2| < cz^ for z larger than some positive constant zq. (2.8) 
Now our main result on the existence of finite energy weak solutions reads as follows. 
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Theorem 2.1. Assume that Vl <Z he a hounded domain with a houndary of class C'^^'^ , 
K > 0, and 7 > |. Then for any given T > 0, the initial- boundary value problem (jl.ip - (|1.2p 
has a finite energy weak solution (p, u, H) on n x {0,T). 

Remark 2.2. The fluid density p as well as the momentum pu should be recognized in the 
sense of instantaneous values (cf. Definition 2.1 in [5]) for any time t e [0, T]. 




As a direct application of Theorem 12.11 we have the following result on the large-time 
behavior of solutions to the problem (|l.ip - p.2p : 

Theorem 2.2. Assume that (p, u, H) is the finite energy weak solution to (|l.ip - ()1.2p ob- 
tained in Theorem \2.1l then there exist a stationary state of density p^ which is a positive 
constant, a stationary state of velocity = 0, and a stationary state of magnetic field 
Hs = such that, as t oo, 

Ps strongly in L'^{Q); 

Us = strongly in £^(ri); (2.9) 
Us = strongly in L'^{fl). 

The proof of Theorem 12. II is based on the following approximation problem: 
' pt + div(pu) = eAp, 

{pu)t + div (pu (K) u) + a\/pt + (5Vp^ + eVu • Vp 

= (V X H) X H + pAu + (A + p)Vdivu, 
^Ht - V X (u X H) = -V X (j/V X H), divH = 0, 

with the initial-boundary conditions which will be specified in Section 4, where /? > is a 
constant to be determined later, and e > 0,S > 0. Taking e — > and (S ^ in ()2.10p will 
give the solution of p.ip in Theorem 1 2. II We remark that the nonlinear term (V x H) x H 
can be dealt with by the idea arising in incompressible Navier-Stokes equations, but there 
are no estimates good enough to control possible oscillations of the density p. In order to 
overcome this difficulty, we adopt the method in Lions [21] and Feireisl [9] which is based 
on the celebrated weak continuity of the effective viscous flux p — (A + 2p)divu. More 
specifically, it can be shown that 

{apZ - (A + 2p)divu„)^(Pn) -> {ap^ - (A + 2p)divu)b{pj 

weakly in L^{n x (0, T)), where p„ and u„ are a suitable sequence of approximate solutions, 
and the symbol F{v) stands for a weak limit of {-F'(wn)}^i- 

Remark 2.3. Similarly to 0, our approach also works for the general barotropic flow: 

p(p) = ap^ + z(p), hm [0,(X3), 

where a > 0,7 > |. Moreover, we also can extend our results to the initial-boundary 
value problem for ()l.ip in an exterior domain by using the method of invading domain (cf. 
Section 7.11 in [24]) and the following special type of Orlicz spaces L^{n) (see Appendix 
A in [21]): 

LP{n) = {fe Ll,{n) /|{|/|<„} e L^in), and /|{|/|>,} e LP{n), for some ry > 0} . 
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3. The Solvability of The Magnetic Field 

In order to prove the existence of solutions to (|2.10p by Facdo-Galcrkin method, we need 
to show that the foUowing system can be uniquely solved in terms of u: 



{Ht - V X (u X H) = -V X (j/V X H), 
divH = 0, (3.1) 
H(a;,0)=Ho, H|ao = 0. 



In fact, we have the following properties: 
Lemma 3.1. Let u e C([0, T]; CQ(ri; M"^)). Then there exists at most one function 

\ieL\[o,T]■,Hl{s^))c^L°-{%T]■,L\n)) 

which solves p.ip in the weak sense on fl x (0, T), and satisfies boundary and initial 
conditions in the sense of traces. 



Proof. Let Hi, H2 be two solutions of (|3.1|1 with the same data. Then we have 

(Hi - H2)t - V X (u X (Hi - H2)) = -V X {uV X (Hi - H2)). (3.2) 



Multiplying ()3.2|) by Hi — H2, integrating over Q, and using the Cauchy-Schwarz inequality, 
we obtain 

~ f |Hi-H2pd2; + i. / |Vx (Hi-H2)pdx 

= / (u X (Hi - H2)) • (V X (Hi - H2))dx 
Jn 

< ^ / |V X (Hi - H2)pdx + ^ / |u X (Hi - H2)pdx 

V f . , , ,0 ,,,,,,/", ,T 

— — — \^A^ I ri \ I xj xj ^, 



< - / |V X (Hi - H2)|Mx + G[v, II u Hoc) / |Hi - H2rdx. 



This implies 



Ydt 



(3.3) 



■ / |Hi - H2pda; + ^ / |V x (Hi - Yi.^)\^Ax 
Jn ^ Jn 

<C(zy,|| ulloo) / |Hi-H2pda;. 
Jn 

Then, Lemma [3T] follows directly from Gronwall's inequality and (|3.3p . □ 



Lemma 3.2. Let 17 C K'^ 6e a hounded domain of class C'^'^'^, k > 0. Assume that 
u G C{[0,T];Cq{^1;R.^)) is a given velocity field. Then the solution operator 

u ^ H[u] 

assigns to u £ C([0, T]; Cq (fJ; R'^)) the unique solution H of (|3.ip . Moreover, the solution 
operator u t-* H[u] maps bounded sets in C([0, T]; Cg (11; R'^)) into bounded subsets of 
Y := L2([0,T];iJi(17)) ni°°([0,r];i2(J7)), and the mapping 

u e C([0,T];Co^(17; R^)) ^HeY 

is continuous on any bounded subsets of C{[0,T];CQ{fl 



Proof. The uniqueness of the solution to ()3.1|) is a consequence of Lemma 13.11 Noticing 
that 

V X (V X H) = V(divH) - AH, 
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then (|3.ip becomes 

fHt V X (u X H) = i/AH, 

I divH = 0, (3.4) 
[h(2;,0)=Ho, H|ao=0, 
which is a hnear parabohc-type problem in H, so the existence of solution can be obtained 
by the standard Faedo-Galerkin methods. And from (|3.3p . we can conclude that the solu- 
tion operator u i— > H[u] maps bounded sets in C([0, T]; Cq (fi; M"^)) into bounded subsets 
of the set Y = L^{[0,T]; H^{n)) n L^{[0,T]; L^{n)). 

Our next step is to show the solution operator is continuous from any bounded subset 
of C([0,T];C2(f];R3) to L^{[0,T];H^n))nL°°{[0,T];L^{n)). To this end, let {u„}^i be 
a bounded sequence in C{[0,T];C^{n;R^), i.e., {u„}^^i C B{0,K) C C([0, T]; C^(17; R^) 
for some K > 0, and 

u„ ^ u in C([0,T];C^(fI;M3)), as n ^ oo. 
Then, we have, denoting H[u] by Hu, and H[u„] by H„, 

~f |H„-Hu|2da- + z. / |Vx (H„-Hu)|2dx 
^ Jn Jn 

= / (u„ X H„ - u X Hu)) • (V X (H„ - Hu))da; 

< f ((u„ - u) X H„ + u X (H„ - Hu)) • (V X (H„ - ll^))dx 

<|| u„ — u||oc-||H„j|y + if ||H„ — Hu||l2(q) II H„ - Hu||_H-i(n) 

< C2||u„ - u||«, + c||H„ - H„||2,(,,) + ^||H„ - Hu||^,(^), 

where, we used the fact that H„[u„] is bounded in Y, says, by C = C{K, v). This implies 
that 

l|/jH..-H„Pd.+ ^/jVx(H,.-H„)Pd. ^^^^ 

< C2||U„ - Ujloo + C||H„ - Hu||i2(o) 

Integrating (j3.5p over t G (0, T), and then taking the upper limit over n on the both sides, 
we get, noting that u„ ^ u in C([0, T]; C^(II; R^)), 



1 limsup / |H„ - Hupd.T + limsup / /" |V x (H„ - Hu)pda;ds 

2 n Jn 2 „ Jo Jn 

< climsup / ||H„ - Hu||i2(f2)ds < c / limsup ||H„ - Hu||^2(f2)ds, 

n Jq Jq n 



(3.6) 



thus, from p.6p . using Gronwall's inequality and the same initial value for H„ and Hu, we 
get 

limsup/ |H„ - Hup(i)da: = 0. 
n Jn 

This yields, from (|3.6p again, 

limsup f /" |V X (H„ - Hu)pdxds = 0. 
n Jo Jn 

Therefore, we obtain 

H„ ^ Hu in Y. 

This completes the proof of the continuity of the solution operator. □ 
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4. The Faedo-Galerkin Approximation Scheme 

In this section, we establish the existence of solutions to (|2.10p following the approach 
in [TOj with the extra efforts to overcome the difficulty arising from the magnetic field. Let 

Xn = span{ryj}"^i 

be the finite-dimensional space endowed with the Hilbert space structure, where the 
functions i]j G M'^), j = 1, 2, form a dense subset in, says, Cg (fl; M"^). Through this 
paper, we use T) to denote C^, and V for the sense of distributions. The approximate 
velocity field u„ e C([0,T];X„) satisfies the system of integral equations: 

/ pvLn{x,t) ■ r]Ax - / mo,5-?7da; 

= / (^Au„ - div(pu„ (g) u„) + V ((A + ^)divu„ - ap''' - Jp'') (^-l) 
Jo Jo. ^ 

- eVp ■ Vu„ + (V X H) X h) • dxdr, 



for any t G [0,T] and any ry G X„, where e, (5, and (3 arc fixed positive parameters. The 
density p„ — p[^n] is determined uniquely as the solution of the Neumann initial-boundary 
value problem (cf. Lemma 2.2 in [10]): 

pt -f div(pu„) = eAp, 

Vp • n\sn = 0, (4.2) 
^p(a;,0) = Pa,s{x); 

and the magnetic field H„ = H[u„] as a solution to the system (|3.ip . The initial data po,* 
is a smooth function in (f2) satisfying the homogeneous Neumann boundary conditions 
Vpo> • n|gn = 0, and 

< <5 < po,5 < <^"^^, (4.3) 
Po^5 ^ Po in L^{n), |{po,5 < Poll ^ for (5 ^ 0; (4.4) 

moreover, we set 

mo,(.) = Ir^'^' I - (4-5) 
0, if po,5 < Pq{x). 



Due to Lemma 13.11 and Lemma 13.21 the problem (|4.1[) , (|4.2[) , and (|3.ip can be solved 
locally in time by means of the Schauder fixed-point technique; see for example Section 
7.2 of Feircisl [9]. As in [lOl [9] the role of the "artificial pressure" term 5p^ in (|4.ip is to 
provide additional estimates on the approximate densities in order to facilitate the limit 
passage e — > (cf. Chapter 7 in [5]). To this end, one has to take (3 large enough, says, 
/3 > 8, and to re-parametrize the initial distribution of the approximate densities so that 

„/3 



pIj^ dx ^ as (5 0. (4.6) 

To obtain uniform bounds on u„, we derive an energy equality similar to (|2.5p as follows. 
Taking 'q{x) ~ u„(x,i) with fixed t in (|4.ip and repeating the procedure for a priori 
estimates in Section 2, we deduce a "kinetic energy equality": 

+ / {p\Du.n\'' + (A + p)(divu„)2 + i.(v X H„) • (V X H„)) dx (4.7) 
Jn 

a7Pr' + '^/3pr')|Vp„pda: = 0. 
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The uniform estimates obtained from (|4.7p furnish the possibiUty to repeat the above fixed 
point argument to extend the local solution u„ to the whole time interval [0,T]. Then, by 
the solvability of equations (|4.2p and (|3.ip . we obtain the functions {p„, H„} on the whole 
time interval [0, T]. 

The next step in the proof of Theorem 12.11 consists of passing to the limit as n ^ oo 
in the sequence of approximate solutions {(0„,u„,H„} obtained above. Wc first observe 
that the terms related to u„ and p„ can be treated similarly to Section 7.3.6 in [3], due to 
the energy equality (j4.7p . It remains to show the convergence of the sequence of solutions 
{H„}^i. From (|tj)) . we conclude 

H„ is bounded in {[0,T]; L^{n)) f] L^{[0,T]; H^{n)). 

This implies that, by the compactness of H^{il) ^ and selecting a subsequence if 

necessary, there exists a function H e L'^{[0,T]; L'^{n)) n L'^{[0,T]; H'^{n)) with divH = 
such that H„(-,i) H(-,t) in L^{n) for a.e. t e [0,r]. Thus, 

(V X H„) X H„ ^ (V X H) X H in V'{n x (0,T)). 

Similarly, we have 

V X (u„ X H„) ^ V X (u X H) in V'{n x (0,T)), 

where 

u„ ^ u weakly in L'^{[0,T]; H^(n)). 
Therefore, (|4.ip and (|3.ip holds at least in the sense of distribution. Moreover, by the 
uniform estimates on u, H and the third equation in p.ip . we know that the map 

t ^ 'H.n{x, t)Lp{x) dx for any (p £ 
Jn 

is equi-continuous on [0,T]. By the Ascoli-Arzela Theorem, we know that 

H(a;, t)ip{x) dx 

is continuous for any ip G V{n). Thus, H satisfy the initial condition in (|3.ip in the sense 
of distribution. 

Now, we are ready to summarize an existence result for problem (|4.ip . (|4.2p . and (|3.ip 
(cf. Proposition 7.5 in 0). 

Proposition 4.1. Assume that C is a bounded domain of the class C^"*"", k > 0. 
Let e > 0, (5 > 0, and (3 > max{4,7} be fixed. Then for any given T > 0, problem (j4.ip . 
()4.2p . and p.ip admits at least one solution p, u, H in the following sense: 
(1) The density p is a non-negative function such that 

p e U{[o,T]-w^^^{n)), dtp e L^((o,r) x n), 

for some r > 1, the velocity u belongs to the class L'^{[0,T]; Hq{Q,)), equation (j4.2p 
holds a.e on n x (0, T), and the boundary condition as well as the initial data 
condition on p are satisfied in the sense of traces. Moreover, the total mass is 
conserved, especially, 

/ p{x,t)dx^ / po,s{x)dx, 
Jn Jn 

for all t G [0,T]; and the following estimates hold: 
5 [ [ p'^'+Mxdt < C(e), 



' / / jVppdxdi < C with C independent of e. 
Jo Jn 
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(2) All quantities appearing in equation (|4.ip are locally integrable, and the equation 
is satisfied in T)' {p, x (0,T)). Moreover, we have 



puGC([0,r];L£^,(f^)), 

and pM satisfies the initial condition. 

(3) All terms in p.ip are locally integrable on O x (0, T). The magnetic field H satisfies 
equation p.ip and the initial data in the sense of distribution. divH = also holds 
in the sense of distribution. 

(4) The energy inequality 

\pu' + -^p^ + + i|Hp) (X, t) dx 

[ {p\Du\^ + A(divu)2 + i/| V X Hp) {x, t) dx 
Jn 

a^p^-^ + 5(3p'^-^) |V/9|^(x,t) dx 
holds a.e t £ [0, T] . 

In the next section, we will complete the proof of Theorem 12.11 by taking vanishing 
artificial viscosity and vanishing artificial pressure. 



(4.8) 



5. The Convergence of the Approximate Solution Sequence 

Now we have the approximate solutions {pe.s, ^e,s^ .5} obtained in Section 4. To prove 
Theorem 12.11 we need to take the limits as the artificial viscosity coefficient e ^ and as 
the artificial pressure coefhcient ^ — > 0. 

First, following Chapter 7 in [9] (see also [10|), we can pass to the limit as e ^ to 
obtain the following result: 

Proposition 5.1. Assume <zM? is a bounded domain of class C^^'^, k > 0. Let S > 0, 

and 

P > max{4, -} 

27 — 3 

be fixed. Then, for given initial data pojUio o.s i'^ (j4.3p - (|4.5p and Hq as in (II. 2p . there 
exists a finite energy weak solution p, u, H of the problem: 

pt + div(pu) = 0, 

{pu)t + div (pu (8) u) + V(ap'^ + (5p^) = (V x H) x H + ^jAu + (A + ^)u, (5.1) 
Ht - V X (u X H) = -V X (i/V X H), divH = 0, 

satisfying the initial boundary conditions (|4.3p - (|4.6p and (|1.2p . Moreover, p G L^^^{Vl x 
(0,T)) and the continuity equation in (j5.ip holds in the sense of renormalized solutions. 
Furthermore, p, u, H satisfy the following uniform estimates: 

sup llp(OlllTm) ^ '^■^'5[po,mo,Ho], (5.2) 

te[o,T] ^ ' 

S sup \\p{t)\\l„,^^. < cEs[po, mo, Ho], (5.3) 

tG[0,T] 

sup \\y^u{t)\\l2m)<cEs[po,nia,B.o], (5.4) 

te[o,T] 

I|u||l2([o,t];h1(^^)) - c£^<5bo,mo,Ho], (5.5) 



GLOBAL SOLUTIONS OF MAGNETOHYDRODYNAMIC EQUATIONS 11 

sup ||H(t)||i.(f,) <cS4po,mo,Ho], (5.6) 
te[o,T] 

l|H||L2([o^T];Hi(f2)) < cEs[po,nio,lio], (5.7) 
where the constant c is independent of 5 > and 

Observing that the conditions (|4.4p - (|4.6p imply that the term Es[po,mo,ilQ] appearing 
in (|5.2p - (|5.7p can be indeed majored to be a constant E[po,m.o,tla] which is also inde- 
pendent of the choice of S. This gives us a list of a priori estimates on p, u, H which are 
independent of S. We omit the proof of Proposition l5.ll and concentrate our attention on 
passing to the limit in the artificial pressure term to establish the weak sequential stability 
property for the approximate solutions obtained in Proposition 15. II as S ^ 0. 



5.1. On the integrability of the density. We first derive an estimate of the density 
ps uniform in (5 > to make possible passing to the limit in the term Spg as (5 ^ 0. The 
technique is similar to that in [lO] . 

Noting that the function b{p) = ln(l + p) satisfies the condition (|2.8p . and ps, us, Hs are 
the solution to (j5.ip in the sense of renormalized solutions, we have 

(ln(l + ps))t + div (In {1 + ps) us) + ( In (1 + ps)] divus = 0. (5.8) 

V 1 + /3<5 / 

Now we introduce an auxiliary operator 

B: l^feL'^m-. lj = 0^^[W^'P{n)r 

which is a bounded linear operator, i.e., 

\\^lf]\\w^-''{n) ^ (^ip)\\f\\Lp{n) for any 1 < p < 00; (5.9) 
and the function W = B[f] E R"^ solves the problem 

divW = / in Q, W\an = 0. (5.10) 
Moreover, if f can be written in the form / = divg for some g G L"^ , g ■ n|ao = 0, then 

\\B[.f]\\Lnn) < c{r)\\g\\Lri,n) (5.11) 

for arbitrary 1 < r < 00. 
Define the functions: 



(fii = i'{t)Bi 



\n{l + ps) - f lii{l + ps)dx 



'i/j e v{o,T), i = 1,2,3, 



where ln(l + ps) dx ~ J^^ ln(l + ps)dx is the average of ln(l + ps) over Q. By virtue 
of ([Q]) and jEH), we get 

ln(l + ps) e C{[0,T];LP{Q)) for any finite p > 1. 

Therefore, from (|5.9p . we have 

ip^ e C{[0,T];W^'P{n)) for any finite p > 1. 

In particular, tpi Cz C{Q x [0,T]) by the Sobolev embedding theorem. Consequently, ipi can 
be used as test functions for the momentum balance equation in (|5.ip . After a little bit 
lengthy but straightforward computation, we obtain: 

/ / ^i5p^,+ap])Hl + ps)dxdt = J2lj: (5-12) 
JoJn -^-^ 
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where 



h = 



tp / {SPs + (^Ps) '^^'f ln(l + Pi) dxdt, 
Jn Jn 



/2 =(A + ^) / / -01n(l + P5)divu5 dxdi, 
JoJn 



h 
h 
h 



\n{l + ps) - / \n{l + ps)Ax 



Axdt, 



ln{l + p5) - f ln{l + ps)dx 



dxdt, 



i^psUgBi 



(ln(l+p.)-^)div 
In(l + P5) 



PS 



l + ps 



divui dx 



dxdt. 



?/>/9i-W5-B,;[div(ln(l + ps)us] dxdt, 



10 Jn 

r-T I 



h = 



V^(V X Hi) X H,- • B, 



Jn 



\n{l + PS ps))dx 



Now, we can estimate the integrals Ii — /y as foUows. 

(1) First, we see that Ii is bounded uniformly in 5, from 
property: 

ln(l + i) 
hm ^ ' = 0. 

(2) As for the second term, we also have 

1-^2 1 < / / \^An{l + ps)diY\is\ dxdt < c, 

by the Holder inequality, (|5.5p . (|5.2p . and the following property: 

ln2(l+t) 



dxdt. 



, and the following 



lim 



0, 



where and throughout the rest of the paper, c > denotes a generic constant. 

(3) Similarly, for the third term, we have 

I/3I < / / i'tPsu^B, 
JoJn 

Here, we have used (j5.4p . (jS.Sp . and the embedding W^-'^'^'(ri) > for p > 3, since 

ln(l + ps) - /j^ ln(l + ps) dx S Lp{Q) for any 1 < p < 00. 

(4) Similarly to (3), we have 



ln(l + - / ln(l + p5)da; 
Jn 



dxdt < 



Jn 



ipdxjUsdx^Bi 



ln(l +ps)-f ln(l + ps) dx 



n 



and, by (|5.2p . (|5.5p . and Holder inequality, we have 



oJn 



IppSUsU^da^^B, 



ln{l +ps) - f ln{l+ps)dx 



n 



dxdt 



dxdt 



< c. 



Here, we used the restriction 7 > |. Therefore, wc obtain 

1/4 1 <c. 
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(5) Next, by Holder inequality and (|5.9p . wc have, 

[■(/'IIIiOaII lk/]OAUAll r2^0lll-B,-MII '2-1 



h\<c \i;\\\ps\\l-,,^A\y/p^U5\\L2(n)\\B^[w]\\ 2. dt < c, 

In ^ ' L { 



since 



w -.^ (\n{l + ps) - — ^ ) divu5 - / ( ln(l + pg) - ) divu^ dx £ L'-(r!), 

V 1 + P5 / Jn\ ^ + PsJ 



for some 1 < ?■ < 2, and here we have used the estimates (|5.2p . (|5.4p . 
(6) Similarly to (5), using (|5.2p . (|5.4p . and we have 

141 < / / \i'PsulB,[dw{ln{l + ps)us)]\dxdt<c. 



./fJ 



Here, we have also used the property (|5.1ip . 

(7) Finally, using Holder inequality again, we have 



41 < c / IV^IIIV X H5|li2(n)||H5|U2(n)dt < c. 
Jo 



Here we used the result Lpi e C{[0,T] x n), ((O)) . and ((5J)) . 
Consequently, we have proved the following result: 

Lemma 5.1. The solutions ps of system (|5.1|1 a^so satisfies the following estimate 



I [ ij{Sp^^ + ap])ln{l + ps)dxdt<c, 
JoJn 



where the constant c is independent of S > 0. 
Remark 5.1. Lemma |5. II vields 



(Jpf + ap])ln(l + ps) dxdt < c. 



Jn 



where c does not depend on S. Using the similar method to Lemma 4.1 in [10| . it can be 
shown (cf. [m [TUl [5] ) that the optimal estimate for the density ps is the following: 



n {6p^,+apJ)pUxdt< 
Jo Jn 



where the constant c is independent of 5 > 0, and > is a constant. But as shown later, 
our estimate in Lemma |5. II is enough for our purpose. 



Define the set 

= {(x, t) e (0, T) X r2 : ps{x, t)<k}, k > 0, 6e (0, 1). 
From (|5.2p . there exists a constant s G (0, oo) such that, for all 5 e (0, 1) and fc > 0, 

meas{r2 x (0, T) - jf} < 



k' 

We have the following estimate: 



Spg dxdt 



< Sp'^g dxdt + Sp^g dxdt 
JJji JJnx(o,T]-Ji 

< Tdkf^ mca.s{n} + S J J Xox(o,T)-,/f Pi dxdt. 
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Then, by the Holder inequahty in OrUcz spaces (cf. [T]) and Lemma l5.ll we obtain 



< 5 (^7V"i {^^ ^ max{l, 2(1 + )hi(l + p^) dxdt} 



< N 



max{l,(41n2)rmeas{17} +4/3 /" [ p^ln(l + pa) dxdt} 

Jo Jsin{p5>l} 

max{(5, (41n2)(5rmeas{ 17} + 4(5/3 /" [ ln(l + pa) dxdt}, 



(5.14) 

where -Ljv/(0), and Ln[Q) are two Orhcz Spaces generated by two complementary N- 
functions 

M{s) = (l + s)ln(l + s)-s, 
N(s) ^ e — s — 1, 



respectively. Due to Lemma |5. 11 we know that, if 5 < 1, 

r - . . < c, 

/o Jn 

for some c > which is independent of S. Combining (|5.13p with (|5.14p . we obtain the 
estimate 



max{<S, (41n2)(5rmcas{r2} + 45/3 [ [ pf ln(l + pa) dxdt} 

Jo Jn 



Spg dxdt 



< Tdk'^ineas{n} + c /V"^ ' ' ' 



10 Jn 

where c does not depend on S and k. Consequently 

lim sup 

a^o 



Jn 



Spg dxdt 



<c[N-'(-)) \ (5.16) 



The right-hand side of (|5.16p tends to zero as fc — > oo. Thus, we have 

lim / / (5pfda;dt = 0, 

which yields 

5p^ ^ in X>'(r2 X (0,r)). (5.17) 

5.2. Passing to the limit. The uniform estimates on p in Lemma [5.11 and Proposition 
15.11 implv. as (5 ^ 0, 

pa ^ p in C([0,T];L;^^^,(^^)), (5-18) 

ua -> u weakly in L^{[0,T]; H^{n)), (5.19) 

and 

Ha weakly* in i2([o,T];i/i(r!))nX°°([0,T];L2(f])), 

divH = 0inI?'(Ox (0,T)); ^ ' -* 
and, from Lemma 1 5. II and Proposition 2.1 in [9], we have, as 5 — > 0, 

Ps-^W weakly in L^([0, T]-L^{n)), (5.21) 

subject to a subsequence. 

By (|5T9l) . ([QO]) and the compactness of Hl{Vi) ^ £^(0), we obtain, 

V X (ua X Ha) ^ V x (u x H) in V'{n x (0, T)), (5.22) 
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and 

(V X Hi) X H5 ^ (V X H) X H in V'{n x (0, T)), (5.23) 

as (5 ^ 0. On the other hand, by virtue of the momentum balance in (jS.ip and estimates 
(|5?2|l -([5J l) . we have, as (5 ^ 0, 



psus ^ pu in C([0,T];L;:^,(f})). (5.24) 
Similarly, we have, as (5 ^ 0, 

H,^HinC([0,T];LLafc(^^))- 

Thus, the limits p, pu, H satisfy the initial conditions of (|1.2[) in the sense of distribution. 

Since 7 > |, ()5.24p and (j5.19|l combined with the compactness of iJ^(ri) ^ i^(ri) 
imply, as (5 — > 0, 

psus (E)us^ pu(E)u in X>'(17 x (0, T)). 
Consequently, letting (5 ^ in (|5.ip and making use of (|5.18p - (|5.24p . (p, u, H) satisfies 

dtp + div(pu) = 0, (5.25) 
dtipu) + div(pu (g) u) - ^jAu - (A + /i)Vdivu + aVp^ = (V x H) x H, (5.26) 
Ht - V X (u X H) = -V X (j/V X H), divH = 0, (5.27) 

in 'D'{n X (0, T)). Therefore the only thing left to complete the proof of Thcorcm l2.1l is to 
show the strong convergence of ps in or, equivalently, pT = p'^ . 

Since ps, us is a rcnormalizcd solution of the continuity equation ()5.ip in I?'(K'^ x (0, T)), 
we have 

Tk{ps)t + dW{Tk{ps)us) + {U{ps)ps - Tk{ps))dWus = in V'{R^ x (0,T)), (5.28) 
where Tk is the cut-off functions defined as follows: 



Tfc(z) = fcT(^-j for z G i?, fc = 1,2,. 



and T £ C°°{R) is concave and is chosen such that 

r(z) = 

Passing to the limit for 5 0+, we obtain 



z, z <1, 
2, z>3. 



dtTkip) + div(Tfc(p)u) + {T;^{p)p - rfc(p))divu = O in 2?'((0, T) x M^)^ 

where 

{Tl{ps)ps ~ rfc(p5))divu5 ^ (T^(p)p-Tfc(p))divu weakly in h'iSl x (0,T))), 

and 

Tk{ps) ^ TMp) in C([0,r];iLafc(f^)) for alll < p < 00. 

5.3. The effective viscous flux. In this section, we discuss the effective viscous flux 
p(p) — (A + 2/i)divu. Similarly to [HKTUIIS], we prove the following auxiliary result: 

Lemma 5.2. Let ps, us he the sequence of approximation solutions obtained in Proposition 
dSU). Then, 

lim / ijj 4>{ap] - {X + 2p)diyus)T,,{ps)dxdt 



/ 0(apT - (A + 2p)divu)Tk{p) dxdt, 
Jn 

for any tp e V{0,T) and (j) S V{n). 
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Proof. As in [lOl [9] , we consider the operators 

AH = ^-^[do^M, « = 1,2,3 

where A^^ stands for the inverse of the Laplace operator on M.^. To be more specific, Ai 
can be expressed by their Fourier symbol 



, i = 1,2,3, 



with the following properties (see [lO]): 



\\Av\\w^--(n) < c(s, r2)||u||L»(fl3), 1 < s < oo 

< c{q,s,n)\\v\\Ls(^ji3-j, q finite, provided - > - - 

q s 6 

\\Av\\L^(^n) < c(s,f^)||v||La(fl3), if s > 3 



Next, we use the quantities 



^,{t,x) ^ xP{t)^ix)A[Tkips)], ^peV{o,T), <f>ev{n), i = 1,2,3, 



as the test functions for the momentum balance equation in (j5.1|l to obtain. 



i'(f>{ap] + 5p1 - (A + 2^)divu5)Tfc(p5) dxdt 



oJn 



-092;, 0( (A + p)dWus - apj - Spf)At[Tk{ps)] dxdt 



+ P f [ V'(5.,^9,,4A:[Tfc(p5)] -49x,09,^A[rfc(P5)] +MR.</'rfc(p5)) dxdt 

Jo Jn 

- [ [ cjypsul {dti^A[Tk{ps)] + M.[{Tk{p5) - U{ps)ps)divus]) dxdt (5.29) 
Jo Jn 



10 Jn 

cT I 



'ippsululdx^(j)At [Tkifs)] dxdt 



ipug (Tk{ps)Ut,j[psul] - (fipsu^gTlt.jlTkips)]^ dxdt 

[ [ V</'(V X H,-) X Hi- • A[n{ps)] dxdt, 
Jo Jn 



10 Jn 



where the operators TZij — dxjAi[v] and the summation convention is used to simplify 
notations. 

Analogously, we can repeat the above arguments for equation (j5.26p and the test func- 
tions 



ip,{t,x) ^ iP{t)(t>{x)A^[Tk{p)], z = 1,2,3, 
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to obtain 

n'il;(l>{ayi' - (A + 2/x)divu)Tfe(p) dxdt 

= 11 V'a..0((A + M)divu - a-^)A^[Tk{p)] dxdt 
JoJn 

+ ^ ^ ^ (a,^ A [T\{p)] - u'a,^ [Jlip)] + u'd^^ mip)) ^xdt 

~ LL '^^^ ^ "^"^^ [(Tfc(p)-r^(p)p)divu]) dxdi (5.30) 

- /" / VP"'"'9,,0A[TUp)]da;dt 

^ i) X "^"^ (j^'^''^ ^^P""'^ ~ (|}pu'n^J [Tfc(p)]) dxdi 

- /" / V0(V X H) X H • yt[7Mp)] d.Tdt. 

Similarly to [TOl E] , it can be shown that all the terms on the right-hand side of ()5.29|) 
converge to their counterparts in (|5.30p . Indeed, with the relations (|5.18p - (|5.24p and the 
Sobolev embedding theorem in mind, it is easy to see that it is enough to show 

J J ^ug (Tk{ps)TZij[(t)psul] - (jipsu^TlijiTkips)]^ dxdt 

^ iX"^"' (tMpI-R^j ['/'P"-''] - (ppu^n^Anip)]) dxdt, 

because the properties of Ai and the weak convergence of u in i^([0, T]; iJ^(fi)) imply 
A^iTkips)) ^ A(TMp)) in C((0,T) x n), 

Tl,ATk{ps)) n^^j{Tk{p)) weakly in LP{[Q,T] x il) for aU 1< p < oo, 

and 

A[{Tk{ps) ~ TUp)p)divus] ^ A^[{n{p)-Tl^{p)p)divu] weakly in L\[0,T]; H\n)). 
From Lemma 3.4 in |10] . we have 

Tk{ps)'Tii,][(j)psul] ~ 4'Psul'R't.][Tk{ps)] 

-^Tk{p)Ri,j[cj)pu^]~ (j^pu^TZi^Tkip)] weakly in i''(r2), i, j = 1, 2, 3, 
for some r > 1 . Hence, we complete the proof of Lemma 15.21 □ 

5.4. The amplitude of oscillations. The main result of this subsection reads as follows, 
and is essentially taken from [1^ (cf. Lemma 4.3 in |10]): 

Lemma 5.3. There exists a constant c independent of k such that 
limsup WTkips) ~ Tk{p)\\L-,+i{{o.T)xn) < c 

5^0+ 
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Proof. By the convexity of functions t p{t), t —Tk{t), one has 
hmsup / / (p]Tk{ps)-7^in{p)) 



T/ 

5^0+ JaJn"" 

hmsup / [ {pj- p'){Tk{p5) - Tk{p)) dxdt 
(5-+0+ Join 



(5.31) 



10 Jn 



{p'f~p''){n{p)~n{p))dxdt 



> hmsup / [pJ- p^){Tk{ps) - Tk{p)) dxdt. 
5^0+ JoJn 

On one hand, we have 

ry ry 

y'' — = / 7s'''^"'"d.s > J (s — zy~^ds = j{y — z)'^ , 



for all y > z > 0, and 

\Tkiy)-mz)\-y <\y-z\\ 

thus, 

(z^ - y-'Kniz) - n{y)) > 7|r,(z) - Tkiy)\->\Tk{z) ~ Tfe(y)| 

= j\n{z)-T,{y)\''+\ 

for all z, y > 0. On the other hand, 



limsup / / [diYUsTkips) - diYuTkip)) dxdt 
5-+0+ JoJn 



limsup / / (Tkifs) - Tfe(p) + Tfc(p) - Tfc(p)) divu^ dxdt 

A^n^ In. la ^ ^ 



5^0+ Join 

< 2 sup i|divu,5||i2((o.T)xn) limsup \\Tk{ps) - Tk{p)\\L'^({o^T)xQ.) {^-^2) 

5 5^0+ 

< dim sup \\Tk{ps) - Tk{p)\\L^((o,T)xQ.) 

(5^0+ 

The relations (|5.3ip . (|5.32p combined with Lemma [5T2l vield the desired conclusion. □ 



5.5. The renormalized solutions. We now use Lemma 15731 to prove the following crucial 
result: 

Lemma 5.4. The limit functions p, u solve (|5.25p in the sense of renormalized solutions, 
i.e., 

dtb{p) + div(6(p)u) + {b'{p)p - fe(p))divu = 0, (5.33) 
holds in I?'(IR^ x (0,T)) for any b G C^{R) satisfying b'{z) = for all z € R large enough, 
say, z > M , where the constant M may depend on b. 

Proof. Regularizing (|5.28|) . one gets 



dtSm[Tk{p)] + div(5™[Tfc(p)]u) + Sm[{T!^ip)p - rfc(p))divu] = rra, (5.34) 

where Sm [v] — Wm *v are the standard smoothing operators and in ( [0, T] ; {M.^ ) ) 

for any fixed k (see Lemma 2.3 in [IS]). Now, we are allowed to multiply (|5.34p by 
b' {Sm[Tkip)])- Letting m — ^ oo, we obtain 



dtb[n{p)] + div(6[Tfc(p)]u) + ib'iTkip))nip) - &(Tfe(p)))divu 

(5.35) 

= 6'(T,(p))[(r^(p)p-Tfc(p))divu] in I?'((0,T) xR3). 
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At this stage, the main idea is to let fc ^ oo in (|5.35|) . We have 



Tk{p) ^ p in LP{n X (0, T)) for any 1 < p < 7, as fc ^ 00, 

since 



\\Tkip) - pI1lp(Ox(o,t)) < hminf \\Tkip5) - P5|1lp(Ox(o,t)), 

and 

mips) P^lli.(ox(o,T)) < 2^^^""IIP^-|II.(nx(o,T)) < cfc^""- 
Thus (|5.35p wiU imply (j5.33p provided we show 



6'(Tfe(p))[(T^(p)p-Tfc(p))divu] ^Oin L\n x {0,T)) as fc 00. 
To this end, let us denote 



Qk,M = {it,x)enx (o,r) I Tk{p) < A/}, 

then 



^' 'b'{Tk{pmmp)p~nip))dwu] 



II {n{p)p-n{p))dwu 



dxAt 

dx&t 



Jo. 

< sup \b'{z)\ 

0<z<M 

< sup \b'{z)\ liminf H(T^(pi-)PA- - Tkip5))divus\\L^Q ^^-^ 

0<^<M d^0+ 

< sup |6'(z)| sup ||uA-||i2([o.T]:_ffi(o)) liminf !|rfc(p5)p5 - '71-(p5)||l2(q 

0<2<A/ 5 A^0+ 

< climinf \\Tl.{ps)ps - Tk{p5)\\ 

o — 

Now, by interpolation, one has 

l|2 

r) 

T + l 

< l|Tfc(P5)P5 -rfc(p5)IU(f2x(0,T))ll^fe(^'5)^'5 -^'^(^<5)lli-; + i(Q,_^,)- 

Similarly to ([5?36l) . we have 

- rfc(p5)||ii(ox(0,T)) < cA:^"''sup||p5||]^^ < cfc^"'^, 



l|rfc(p5)pi-Tfe(p5)||i2(Q^ 



and, using T^(z)z < Tk{z), 
\\\Tk{Ps, 

< I|rfe(p5) - rfc(/?)||i^+i(ox(o,T)) + ||rfc(p)|ji^+i(Q^^ 

< l|T'fc(P<5) - Tk{p)\\L-' + i{nx{0,T)) + \\Tk{p)\\L-' + ^Qk,M) 

+ \\Tkip) - Tk{p)\\Li+^nx{o,T)) 

< \\Tk{ps) - Tfe(p)||L-'+i(nx(o,T)) + Mc{n) 

+ W^kip) - Tk{p)\\Li + i{nx{0,T))- 

From Lemma [5.31 and ()5.38p . we obtain 

\imsuv\\Tiips)ps-Tkips)\\L~,^iiQ, < 4c + 2Mc(r!), 

.5^0+ 



which, together with (|5.37p - ()5.38p . completes the proof of Lemma \5A 
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5.6. Strong convergence of the density. Now, we can complete the proof of Theorem 
12.11 To this end, we introduce a sequence of fimctions Lk G C^{R): 

j zlnz, < z < k, 

'^^^ " \zHk) + zj^ ^ds, z > k. 

Noting that Lk can be written as 

Lk{z)^ Pkz + bk{z), (5.39) 

where bk satisfies the conditions in Lemma 15.41 we can use the fact that ps,us are renor- 
mahzed solutions of (|5.ip to deduce 

dtLkips) + dW{Lk{ps)us) + Tk{ps)dWus = 0. (5.40) 

Similarly, by (|5.25p and Lemma [5T31 we have 

dtLkip) + div(Lfc(p)u) + Tfe(p)divu = 0, (5.41) 

in X''((0,T) X n). By ((OO)) . we can assume, as (5 ^ 0, 



Lkips) Lkip) in C{[0,T];Ll^^,m. 
Taking the difference of (|5.40[) and (|5.4ip and integrating with respect to t, we get 

{Lk{ps) - Lk{p))<t>dx 

t . (5.42) 

/ {{Lk{ps)\is - Lk{p)n) ■ V0 + {Tk{p)diYU - Tk{ps)diYUs)4>) dxdt, 
/o Jn 

for any e 'D{Vl). Passing to the limit for (5 ^ and making use of (|5.42|) . one obtains 



(Lfe(p) - Lk{p))(t)dx 

n 

t 



= / f {Lk{p)-Lk{p))u-Vcl>dxdt (5.43) 
Jo Jn 

+ lim / / {Tk{p)diYU-Tk{p5)diYUs)(t>dxdt, 

for any (j) G 15(51). 

Since the velocity components u', z = 1, 2, 3, belong to i^([0, T]; Wq ' (fl)), one has the 
following (see Theorem 4.2 in P]): 

eL'{[0,T];L\n)). 



distfe, dCl] 



Let us consider a sequence of functions G 'D{Q) which approximate the characteristic 
function of CI such that 

< (pm < 1, 4'm{x) = 1 for all X such that distfa;, dil] > — , 

m (5.44) 

and |V0m(a;)| < 2m for all x ^ il. 

Taking the sequence (f) = (pm as the test functions in (j5.43[) . making use of the boundary 
conditions in (|1.2p . and passing to the limit as m — > cx3, one has 

/ (Uip) - Lk{p)) jj Tk{p)diYndxdt - lim / / rfc(p5)divu5 dxdt. (5.45) 

Jn JoJn ^^^+ JoJn 

We observe that the term Lk (p) — Lk (p) is bounded by (|5.39|) . 

At this stage, the main idea is to let A: — ^ cx) in (|5.45|) . By (|5.2p . we can assume 



PeMPe) pMp) weakly star in L°°([0, T]; L"(fl)) for all 1 < a < 7. 
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Wc also have 



Lk{p) pln(p) in L°°{[0,T];L°'{n)) as fc -> oo for all 1 < a < 7, 
since, by (|5.2p . 

lim r{k) ^ 0, where r{k) := meas{(a;,t) <E fl x {0,T)\ps(x,t) > k}; 

k — ^00 

and because Lk{z) < zlnz, repeating the similar procedure to (j5.14p . we have 



\\Lk{p) - pln(p)!|L-([o.T];L°(n)) 

< sup liminf ||Lfe(p5) - p5ln(p5)||L^([o,T]:L°(n)) 
te[o,T] 5^0+ 



<2g(/s)sup sup max{l, / A'/(p^|lnp5|") dx} 
<5 te[o.T] Jn 

<2g(fc)sup sup niax{l,2 / (1 + p?|lnpa|")ln(l + p^|lnp5|") dx} 
<5 te[o,T] Jn 

<2q{k)sup sup niaxjl, c(a)meas{0} + c(q!) / pg\\nps\°'^^ dx} 

s te[o.T] Jnn{ps>e} 

<2,(fc)sup sup niax{l,c(a)nieasm + c(a,7) /Pjdx} 
s te[o.T] Jq 

< cq{k) ^0, as fc — > 00, 
where the function M is defined in (j5.15p . c is a constant independent of S and 

q{k) \\X[ps>k]h.m < {n~' {;^)) ' 

Similarly, we have 

Lfe(p) ^ pln(p) in L°°([0,r];L"(17)) as fc ^ cx), for alll < a < 7, 
and, by Lemma [5751 



Tk{p) ~*Tk{p) in L"{[0,T];L"{n)) ask ^00, for alll < a < 7 + 1. (5.46) 

Finally, making use of Lemma l5.2l and the monotonicity of the pressure (see ()5.3ip ). we 
obtain the following estimate on the right hand side of (|5.45p : 



/ / rfc(p)divud.Tdi — lim / / Tk{ps)divus dxdt 
Jo Jn '5^0+ Jo Jn 

< I I {n{p)-Tk{p))diYndxdt. 
Jo Jn 



(5.47) 



From (|5.46p and the Sobolev embedding theorem, we see that the right hand side of (|5.47p 
tends to zero as fc ^ 00. Accordingly, one can pass to the limit for /e ^ 00 in (|5.45p to 
conclude 

^ (pHp) - pHp)) (x, t) = 0, for t e [0, T] a.e. (5.48) 
Because of the convexity of the function z zlnz, we have 



pMp) > pln(p), a.e. in n x (0,T), 
which, combining with (j5.48p . implies 



pln{p){t) = pln{p){t), for t G [0, T] a.e. (5.49) 
Theorem 2.11 in [H], combined with (|5.49p . implies 

Pe p, a.e. in X (0, T). 
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From the estimate (|5.2|) on p, together with Proposition 2.1 in [9], again we know, 

ps^ p weakly in L^{n x (0,r)), 

subject to a subsequence. By Theorem 2.10 in [9], we know that for any > 0, there exists 
(T > such that for all (5 > 0, 

/ ps{t, x) dxdt < 77, 
Je 

for any measurable set E C ^ x (0, T) with meas{i?} < a. 

On the other hand, by virtue of Egorov's Theorem, for tr > given above, there exists 
a measurable set C x (0, T) such that 

meas{i?cr} < cr, and ps{x, t) —^ p[x, t) uniformly in H x (0, T) — E^- 

Therefore, we have 

/ / \ps- p\ dxdt 

J jnx{o,T) 

< II \p5-p\dxdt+ 11 \ps-p\dxdt (5.50) 

< 277 + Tmeasjrj} sup \{ps — p){x,t)\, 

{x.t)eE- 

which tends to zero if we first let 5 0+, and then let ?/ — > 0+. The strong convergence 
of the sequence ps in L^(il x (0,T)) follows from (|5.50p . 
The proof of Theorem 12. II is completed. 



6. Large-Time Behavior of Weak Solutions 



Our final goal in this paper is to study the large-time behavior of the finite energy weak 
solutions, whose existence is ensured by Theorem 12. II 
First of all, from Theorem 12. !( we have 



ess supt>o£^(i) + / I {p\Du\^ + {\ + p){d\Ynf + v\'^ xn\^] dxdt <E{fd). (6.1) 
Jo Jn 

Following the idea in [llj . we consider a sequence 

Pm.{x,t) := p{x,t + m); 
Um{x, t) := u(x, t + to); 
Hjn{x,t) := H(.T,i + m), 

for all integer m, and t E (0,1), x E ^l. It is easy to see that (|6.1I) yields uniform bounds 
of 

p,n e L°°([0,l];LT(f])), H,„ e L°°([0,l];L2(f])) 

V^u,„ e L°°i[0,l];L\n)), p^Um e L°°([0,l];i^(r!)), 
which are independent of m. Moreover, we have 



lim f (||Vu„J||.(j,) + ||V X H„ 



2 



dt = 0. 



(6.2) 



Hence, choosing a subsequence if necessary, we can assume that, as to 00, 

Pm{x,t) — > ps weakly in L'^{Q x (0, 1)); 
Ura{x, t) Vis wcakly in ^^([0, l];H^{n)); 
H™(x,0 -> H, weakly in £^([0, 1]; i7o'(f^))- 
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Furthermore, 



/ Psdx < hminf / pm{t) dx < C{Eq). 
Jn Ja 

Therefore, from the Poincarc incquahty and (|6.2p . we know 



This, combined with the compactness of ^ L^, imphes 

Us = 0, a.e in X (0, 1). 

Similarly, we know that 

Hs = 0, a.e in r2 X (0,1). (6.3) 
On the other hand, by Sobolev inequality, Holder inequality, (|6.ip and (|6.2p . we have 



lim / ||/9„|u™p|| jiy + ||/Om|um|||^ 6-, ]dt = 0. 



(6.4) 



Since p, u are solutions to (1.1) in the sense of renormalized solutions, one has, in particular, 

pt + div(pu) = in ^'(f^ X (0,T)). (6.5) 

Then taking test functions Lp{x, t) = ip{t)(f){x) in (|6.5p . where tp{t) G 2?(0, 1), </> G ^{il), we 
have, using integrating by parts, 

[ ( [ p„,(j)dx] il}'{t)dt+ [ [ p^u^V^V dxdt = 0. 
Jo \Jn / Jo Jn 

Letting to — > oo and using (|6.4p . we obtain 

^ ps<l>dx\ tj/{t)dt = 0. 



This implies that p^ must be independent of t, by the arbitrariness of ip. 

Now, following the procedure used in Section 5, or in Lemma 4.1 in [TU], one can obtain 

p'J^^ is bounded in L^{fl x (0, 1)), independently of m > 0, 

for some 9 > 0. Consequently, one has 

p;; -^p^ weakly in l1(J1 X (0,1)). (6.6) 

Therefore, passing to the limit in the momentum balance equation of (jl.ip and using (|6.2p . 
(EH), we get 

Vp^ = in ^'(f^). (6.7) 

Now, we show that the convergence in (|6.6p is indeed strong. To this end, similarly to 
[TT] . we consider 

G(z) = 2;", < a < min i — , 



, 27 ' 61 + 7 J ' 

so that b{z) = G{z'') may be used in (|2.6p . Consider the vector functions 

[G(p7J, 0,0,0] and [p;^„ 0,0,0] 

of the time variable t and the spatial coordinates x. Using (|2.6p and (|6.2p . (|6.4p . we get 

Div[G(p2i), 0, 0, 0] is precompact in W'l'"^' (Q x (0, 1)), (6.8) 

for some gi > 1 small enough. Similarly, making use of the momentum balance equation 
in (fTTT]) . and (IQ)) . we obtain 

Curl[p7„, 0, 0, 0] is precompact in Wi^J'"*" (fi x (0, 1)), (6.9) 
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for some (72 > 1, where 

Div(fo, fi, f2, fa) := (fo)t + S?=i9x.fi, 

and 

Curl(fo,fi,f2,f3) := 9ifj - 9jfi, xq := t, i,j = 0, ...,3. 
Meanwhile, we can assume 



G{p7a) ^ G(p7) weakly in L^- [Q x (0, 1)), (6.10) 

and 

GifDpln ^ weakly in L^{Q x (0, 1)), (6.11) 

with 

1111 
P2 = -, — + — = -<!. 
a pi P2 r 

Using the L^'-version of the celebrated div-curl lemma (see [26]), we deduce that from 

dOjl- ljOTl) 

G{p^)7f=G{p^)p'r. ^ (6.12) 

As G is strictly monotone, (|6.12p implies G{p'^) = G{p'^). Since L~ is uniformly convex, 
this yields strong convergence in (j6.6p . Therefore, we have 

Pm ^ Ps strongly in L'^ [Vl x (0, 1)). 

This, combined with (|6.6p and (|6.7p . gives VpJ = in the sense of distributions, which 
implies that ps is independent of the spatial variables. 

Finally, by the energy inequality, the energy converges to a finite constant as t goes to 
infinity: 

:=Ti^i;(<), 

t^oo 

and, by and 

/•rn+l 

lim 



im / / p|upda; = 0, 
J,n Jn 

/•m+l p 

lim / / IHpdcc^O. 



Thus 



m — -'OQ 



rn J il 

p1 Ax. 



1 9 a „, 1 , 



E^^ lim / / [-pu' + -p^ + -\H.\^] dxdt 



2'^ 7- r 2' 



7-1 

Furthermore, using the continuity equation in (1.1), one easily observe that 

p(x, t) Ps weakly in L^{Cl) as t — > 00. 

Thus, we have 

Eoa = / Ps < liminf / p'^ dx < lim sup / p'^ dx 

Jfi 7 - 1 Jo 7 - 1 t-*oo Jo 7 - 1 

< TSr / f ipu2 + ^-p^ + l\n\A dx = '^E{t) - E^. 
t— ►00 Jjj\z 7— 1 2 / t^oo 



This implies 



lim / — ^p'^ da: = / —^p] dx, 



t- 

and (|2.9p follows since the space L'^ is uniformly convex. 
This completes the proof of Theorem [ 
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